Pseudogap and antiferromagnetic correlations in the Hubbard model 
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Using the dynamical cluster approximation and quantum monte carlo we calculate the single- 
particle spectra of the Hubbard model with next-nearest neighbor hopping t' . In the underdoped 
region, we find that the pseudogap along the zone diagonal in the electron doped systems is due 
to long range antiferromagnetic correlations. The physics in the proximity of (0, n) is dramatically 
influenced by t' and determined by the short range correlations. The effect of t' on the low energy 
ARPES spectra is weak except close to the zone edge. The short range correlations are sufficient 
to yield a pseudogap signal in the magnetic susceptibility, produce a concomitant gap in the single- 
particle spectra near (it, tt/2) but not necessarily at a location in the proximity of Fermi surface. 



Introduction - The normal state phase of high T c su- 
perconductors at low doping, the pseudogap (PG) region, 
is characterized by strong antiferromagnetic (AF) corre- 
lations and a depletion of low energy states detected by 
both one and two-particle measurements^}. Whereas the 
d-wave superconducting phase appears to be universal in 
the cuprates 0,0, the PG region displays different prop- 
erties in the electron and hole doped materials 0, [3. In 
order to further develop a theory for high T c supercon- 
ductivity it is essential to have a better understanding of 
the asymmetry and similarities between the electron and 
the hole doped materials. 

In the hole doped cuprates the antiferromagnetism 
is destroyed quickly upon doping (persisting to « 2% 
doping) 6] and the angle resolved photoemission spec- 
tra (ARPES) show well defined quasiparticles close to 
(-7r/2,7r/2) in the Brillouin zone (BZ) and gap states 
in the proximity of (0,77)0) 0, HI- In the electron 
doped cuprates AF is more robust (persisting to w 15% 
doping) and the ARPES at small doping (w 5%) shows 
sharp quasiparticles at the zone edge and gap states else- 
where in the BZ[j| @- In the Hubbard model, or the 
closely related t-J model, the electron-hole asymmetry 
can be capture d by including a finite next-nearest neigh- 
bor hopping £'|lfj, yT|. In this Letter we employ a re- 
liable technique, the dynamical cluster approximation 
(dca)|H Ea, on relatively large clusters, to investi- 
gate the PG and single-particle spectra at small doping, 
the asymmetry between electron and hole-doped systems, 
and the role of AF correlations on the PG physics. 

We find that in the hole doped systems, the PG emerges 
in the proximity of (0,7r), requires only short range cor- 
relations, and its magnitude and symmetry is strongly 
influenced by t'. In the electron doped systems, the PG 
emerges along the diagonal direction, as a direct con- 
sequence of AF scattering, and requires long range AF 
correlations, but not necessarily long range order. The 
hopping t' enhances the AF correlations in the electron 
doped system and produces this AF gap. With reduced 
temperatures, the short range AF correlations suppress 
the low-energy spin excitations in both electron and hole 



doped systems concomitant with the development of a 
single-particle gap in the proximity of (tt, tt/2), but not 
necessarily with a PG close to the nodal or the antinodal 
points. 

Formalism - The Hubbard Hamiltonian is 

# = E E c U CT +c/E n *T n *f • w 

((H)), <r 1 

Here c^J destroys (creates) an electron with spin a on 
site i and rii a is the corresponding number operator. U 
is the on-site repulsion taken U — W — 8t and t (f) 
is the (next) nearest-neighbors hopping. W is the elec- 
tronic bandwidth. We keep the filling n < 1 and take 
t' = -0.3i (f = 0.3£)[3 to represent the hole (electron) 
doped cuprates. We present results for n — 0.95. 

In the DCA^3| we map the original lattice model onto 
a periodic cluster of size N c — L c x L c embedded in 
a self-consistent host. The correlations up to a range 
£ < L c are treated accurately, while the physics on longer 
length-scales is described at the mean-field level. The re- 
duction to an effective cluster model is achieved by coarse 
graining the BZ into N c cells (see Fig. 3 in Ref.[l3j) and 
approximating the self-energy as a constant within each 
cell, £(k, u>) w £(K,u;), where K denotes the center of 
the cell which k belongs to. 

We solve the cluster problem using quantum Monte 
Carlo (QMC) lla. We use two different 16 site cluster 
geometries [la fl7j . 16A and 16B (see Fig.5 in Ref. 13]), 
which result in different coarse graining of the BZ. Calcu- 
lations on larger clusters below the PG temperature and 
at large coupling (U = W) are not currently possible 
due to the QMC sign problem. The Maximum Entropy 
method^] is employed to calculate the real frequency 
cluster Green's function from which the self-energy is ex- 
tracted. The self-energy is interpolated using a smooth 
spline, and used to calculate the lattice spectrum A(k, uo). 
We find results identical to within error bars at all the 
common points of the coarse-grained BZ of the clusters. 
The result demonstrates that these 16 site clusters cap- 
ture the momentum dependence of the self-energy rather 
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FIG. 2: A{k,Ld) for different k points in the BZ for hole (full 
line) and electron (dashed line) doped cases, at 5% doping 
and T = 0.12t. In (f) the dotted line represents the spectrum 
for t' = 0. 



FIG. 1: (color) 5% doping, T — 0.12f. Zero energy surface 
A{k,0) for a) t' = -0.3t and b) t' = 0.3t. A{k,u>) for k along 
the (0, 0)-(7r, 7r) direction in the BZ for c) t' — — 0.3t and d) 
t' = 0.3i. 

well. We checked the robustness of our results at low 
temperature with calculations on smaller clusters where 
the sign problem is less significant. 

Results - At a temperature T N = 0.19t (0.24t) for the 
hole (electron) doped system the AF correlation length 
reaches the cluster size yielding a divergent AF suscepti- 
bility (not shown). Below TV one can proceed either by 
imposing the full symmetry on the effective medium, i.e. 
by reducing the problem to a cluster embedded in a para- 
magnetic (PM) host, or by allowing the host to develop 
long-range AF order. Both the PM and the AF solutions 
are complementary approximations to the exact solution: 
the first cuts off the AF correlations larger than the clus- 
ter size while the second introduces long range AF order 
via the mean-field character of the host. 

Paramagnetic solution - In Fig. ^ -a and -b we show 
the spectral intensity at zero energy for the hole and elec- 
tron doped systems, respectively. These plots are similar 
to the experimental ARPES data (see Fig. 8 in Rcf. 4] 
and Fig. 3 in Ref. 5]). In both experiment and in our re- 
sults, a region of large intensity can be observed close to 
(7r/2, 7t/2) and very low intensity is observed at the zone 
edge for hole doped systems. For the electron doped sys- 
tems the intensity is maximum at (0,7r). However, the 
experimental data for the electron doped materials show 
gapped states along the diagonal direction Q . In our cal- 
culations, Fig. n]-b, the intensity at (tt/2,tt/2) is similar 
to the one observed for the hole doped case and there is 
no PG along the zone diagonal. In fact A(k, ui) along the 
diagonal direction shows similar features for the hole and 
electron doped cases (Fig. ^-c and -d). Apart from the 



differences at high energy close to the zone center and 
the zone corner which follow the non-interacting disper- 
sion, the low-energy features along the zone diagonal are 
almost identical. 

A comparison of the hole and electron spectra is pre- 
sented in Fig. Inland Fig.|21-a, where A(K, ui) for K in the 
center of the cells which divide the BZ [26j are shown. In 
Figs.|2-a, -b, -c and Fig.|21-a, we find that the single par- 
ticle spectra at low energy for the hole and the electron 
doped cases are surprisingly similar apart from the fea- 
tures close to (0, 7r). In Fig.[21-c we observe a sharp peak 
at (7r/2, 7r/2) in both the hole and electron doped spectra. 
Thus, there is no PG along the diagonal direction 
A particularly interesting feature, shown in Fig. [31 -a, is 
the depletion of the low energy states with lowering T 
in the proximity of (-7T, 7t/2). Unlike the non-interacting 
case, where at (7r,7r/2), there is only spectral weight for- 
ce; > 0, there is now a broad feature with substantial 
weight at negative energies. This is due to AF scattering 
as can be deduced by comparing the main features with 
the (7t/2, 0) spectrum (Fig. |3-a) found at the mirroring 
position with respect to AF zone boundary in the BZ. 
These shadow states develop a gap with decreasing tem- 
perature as shown in Fig.[3]-a where a large temperature 
spectrum (T — 0.22i, dotted line) and a low temperature 
one (T — 0.12i, dashed line) are plotted for the electron 
doped case. In literature the common description^] of 
ARPES along the (it, 0) — (it, ir) cut is that the peak 
at (ir, 0) and u> < characterizing the PG evolves into a 
broad feature which loses intensity when approaching the 
zone corner. Our results indicate that the broad feature 
and gap at (ir, ir/2) are not conditioned by the (0, ir) PG, 
being present in both electron and hole doped systems. 

Differences between the hole and electron doped spec- 
tra are illustrated in Figs. [21 -d through -f. The high 
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FIG. 4: a) Static spin susceptibility versus T. b) A(k,u) 
averaged over the (7r/2,7r/2) cell for the AF and the PM so- 
lutions. 



FIG. 3: 5% doping, hole doped (full line) and electron doped 
(dashed and dotted lines) cases, a) A(k, uj) at k — (tt, tt/2). b) 
Low-energy dynamical spin susceptibility \" (q,u) /u) at q = 
(0, 7r/2). The suppression of the spin excitations starts at the 
same temperature with the depletion of DOS at (tt, 7r/2), see 
dotted lines. c)x"(q,u)/uj at q = (tt/2, tt/2). d) RPA results 
for x"(?,w)/w at q = (tt/2, tt/2). 

energy features at zone center (Fig^-d) and zone corner 
(Fig [21 -e) are strongly influenced by t' . It is interest- 
ing that the position of these features seems to follow 
the non-interacting band structure, the energy difference 
between the non-interacting states at these points being 
about 8|t'|. A fundamental difference between the elec- 
tron and hole doped spectra at (0, n) is shown in Fig. [5} 
f. The hole doped spectra exhibits a strong gap whereas 
the electron doped spectra has an intense peak. It is also 
worth looking at the t' = case, shown in the figure 
with dotted line, where a PG is present but much less 
developed than the one for t' = —0.3t. Notice that even 
for the hole doped case (i.e. t' < 0) the magnitude of t' 
has a strong influence on symmetry with respect to zero 
energy of the density of states, which we believe to be 
important for interpreting tunneling experiments. 

Aside from the depletion of DOS at the chemical po- 
tential, the PG is also associated with the suppression 
of the low-energy spin excitation For both electron 
and hole doped cases we find that the static spin suscep- 
tibility at q = (0, 0) is strongly suppressed at tempera- 
tures below T* w 0.24i, see Fig. @Ja. The momentum 
dependence of the dynamical spin susceptibility exhibits 
a dispersion similar to magnon one in the undoped anti- 
ferromagnets, with zero energy excitations at q = (0, 0) 
and q = (tt, tt), and gapped excitations at other q points 
in the BZ. For instance in Fig. [31 -b we show the imagi- 
nary part of the dynamical spin susceptibility, x"(<7> u)/w 
at q = (0, 7r/2). In the electron doped case the posi- 
tion of the peak is found at larger energy which can be 
interpreted as a larger effective exchange interaction J 
and is consistent with stronger AF. What is interesting 
is that the suppression of the spin excitations and the 



formation of the remnant magnon peaks in x"(<?> w)/cj is 
concomitant with the development of the gap in the sin- 
gle particle spectrum at (tt,tt/2) (see Fig. [3] -a) and not 
necessarily imply a PG at (0,7r), observed only in the 
hole doped case, or elsewhere in the proximity of Fermi 
surface. However, the hole doped PG at (0, tt) is also 
coincident with the appearance of the remnant magnon 
peaks, indicating a short range AF correlations origin. 

The spin excitation spectra are qualitatively similar for 
the electron and hole doped cases, but with a stronger 
suppression at low-energy in the hole doped case even 
though the AF is weaker. The largest difference can 
be observed at q — (tt/2,tt/2) where at T = O.I2i the 
hole doped spectra show a well developed gap whereas 
the electron doped one just starts to form (see Fig. 
-c). This difference is a result of the corresponding dif- 
ferences in the ARPES, as can be concluded from Fig. [3] 
-d. Here the Random Phase Approximation (RPA) using 
the calculated A(k,u>) was employed for the calculation 
of the spin susceptibility. In this approximation the hole 
doped spin susceptibility at (7r/2, 7t/2) is gapped due to 
the gap at (0, tt) in the DOS which suppresses the excita- 
tions between the antinodal and the nodal points, unlike 
the electron doped susceptibility which is peaked at low- 
energy. 

Antiferromagnetic solution - In Fig. 01-b, we compare 
A(k,u) for AF and PM cases. Here, a gap is obtained 
for the AF electron doped case close to (7r/2, it/2) in the 
BZ, in agreement with the experimental findings 5]. This 
gap is an AF gap and requires long range AF correla- 
tions. The short range AF correlations, of the order of 
a few lattice constants, are not sufficient to produce it. 
However, it is possible the PG to also appear in the PM 
state in large enough clusters that allow for long-range 
AF correlations. This conclusion is similar to the weak 
coupling PG mechanism predictions |2lj. even though in 
our case U — W. The hopping t' enhances the antifer- 
romagnetism in the electron doped systems, producing 
the gap. Presumably any other parameters which favor 
the antiferromagnetism will have a similar effect. For ex- 
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ample, the AF solution for the hole doped case produces 
a gap at {tt/2, tt/2) too, though a little smaller due to 
weaker antiferromagnetism. The spectral features away 
from (7r/2, 7r/2) within the AF solutions are not quali- 
tatively different from the ones obtained with the PM 
solution (not shown). We note that the long range AF 
order does not yield a gap at (0, tt) for the electron doped 
case even though this point is on the AF zone boundary. 
The gap in DOS at (tt, tt/2) developed in the PM solu- 
tion is now enhanced by AF order, as well the intensity 
of the shadow states. 

Our conclusions about the nature of the PG in the elec- 
tron doped systems are different from those drawn from 
cluster perturbation theory (CPT)[2l|. For U&W, CPT 
finds that, even when only short range AF correlations 
are considered, the states along the diagonal direction de- 
velop a gap, which persists even at large dopings w 15%. 
Whereas in the experiment, at this doping value, the PG 
shows only at the intersection of the AF zone boun dary 
with the non-interacting Fermi surfaces (hot spots) [22J. 
For agreement with experiment the authors of Ref.|21| 
proposed two different mechanisms for the PG in elec- 
tron doped systems: a strong-coupling (U w W) PG at 
small doping produced by short range correlations and a 
weak-coupling (U < W) PG valid at intermediate doping 
which requires long range AF correlations. In contrast, 
we find no PG along the diagonal direction in the strong 
coupling regime unless long ranged AF correlations are 
considered, implying no need for two different PG mech- 
anisms in the electron doped systems. A plausible reason 
for the discrepancy between the DCA and CPT results 
is the overestimation of AF correlations in the latter ap- 
proach due to finite size effects (lij on small clusters [24j 
and lack of self-consistency j^. 

We find that the inclusion in our calculation of a next- 
ncxt-ncarcst-ncig hbor hopping t" « 0.2*13 will not 
change the conclusions, this term having a rather small 
quantitative effect (though it may provide better agree- 
ment with experimental data). With decreasing doping, 
T* increases but the number of available low-energy un- 
occupied states becomes smaller and therefore the PG 
features are more difficult to be resolved. With increas- 
ing doping, T* decreases and the PG weakens, its fea- 
tures being hardly discernible above 15% doping. In the 
same time the lu < weight in the DOS at (ir, tt/2) is re- 
duced, indicating weaker AF scattering with increasing 
doping. In the AF solution the PG at larger doping will 
be located at the "hot spots", even though the physics 
at those points requires fine k resolution and therefore is 
obtained from interpolation. 

Conclusions - Using DCA we investigated the Hub- 
bard model with next-nearest neighbor hopping t' . We 
find that: (1) The PG along the diagonal direction of the 
BZ in the electron doped systems is an AF gap which 
requires long range AF correlations. (2) The DOS in the 
proximity of (0, tt) is determined by the short range AF 



correlations and it is strongly influenced by t' . For the 
hole (electron) doped systems t' yields a gap (an intense 
peak) at the zone edge. (3) Except in the proximity of 
(0, tt) the influence of t' on the low-energy ARPES is 
very weak. (4) t' has a strong influence on the high en- 
ergy ARPES close to zone center and zone corner. (5) 
The magnitude of t' in the hole doped systems influences 
strongly the symmetry of the PG around the chemical po- 
tential. (6) The short range AF correlations suppress the 
low-energy spin susceptibility, produce remnant magnon 
peaks in the spin excitation spectra in both electron and 
hole doped systems and produce a gap in ARPES around 
(tt, tt/2) but not necessarily in the proximity of the Fermi 
surface. When a (0, tt) PG is present in the ARPES, it 
emerges at the same temperature as the magnon peaks, 
suggesting a common origin of short ranged AF order. 
(7) Even though the antiferromagnetism is stronger in 
the electron doped case, the intense peak in DOS at (0, tt) 
hinders the suppression of low-energy spin excitations. 
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